FLYWHEEL TIMING GENERATION METHOD AND APPARATUS FOR TDMA 
SATELLITE COMMUNICATIONS SYSTEM 

Field of the Invention 

5 This invention relates generally to satellite communications systems, and, more 

pcirticularly, to a flywheel timing generation method as applicable to TDMA satellite system 
design and operation. 

BACKGROUND 

In time-division multiple access (TDMA), terminals take turns using an entire 

10 transmission channel. The terminals transmit according to a frame consisting of a number of 
time slots, each terminal able to use the entire frequency band of the channel during its 
assigned time slot. The time slot is measured from a frame marker which repeats at a fixed 
period, although a time slot can be either fixed or variable in length. The terminals assemble 
packets for transmitting during their assigned time slots. The transmission can be in the form 

15 of bursts having a length of an integer number of slots, the bursts typically consisting of a 
preamble, a xmique word and random symbol data (message portion). 

Delay times, or guard times, are placed between the individual slots to assure that the 
signals from different terminals do not overlap. When the terminals communicate with a 
common terminal or satellite, the terminals need to be synchronized with the satellite at a 

20 tolerance within a fraction of the guard time. The fransmission and reception can be on a 
same channel or on different channels. In a system, multiple terminals communicate with 
each other to synchronize and correct their assigned start times and delay times, to assure that 
the time slot allocation is accurately maintained. 

A receiving terminal is conventionally used to extract a system clock, and 

25 synchronization and timing correction information from a signal transmitted by a satellite. 
The terminal may include a timing controller that generates system timing based on received 
signals from the satellite. The timings are usually based on integer multiples of a number of 
symbols, a symbol being an encoded modulated piece of a larger signal that represents a 
predetermined number of bits of information. Flywheel timing establishes a synchronization 
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time reference for individual traffic terminals. The synchronization involves a time marker 
that the traffic terminal uses to align its transmission, the flywheel timing predicting the 
presence of the synchronization word in received signals based on an anticipated repetition 
frequency and data pattern for the sjTxchronization word. A terminal may conventionally 

5 adjust its system clock when a new synchronization word is received so that the timing of 
transmitted signals is corrected with respect to the satelUte reference clocks. When a flywheel 
circuit receives the synchronization word at the anticipated times, it determines that 
synchronization has been established and sends a reference pulse used to control 
transmissions from individual traffic terminals. The timing confroUer generates a system 

10 clock based on the received signals from the satellite and the reference pulse from the 

flywheel circuit. Either a symbol clock reference or timing correction information is received 
from the satellite. 

TDMA system timing is generated using an adjustment for satellite ephemeris, or 
position information. The range change between the satellite and a terminal is conventionally 

15 predicted using the ephemeris. Then the TDMA timing of the terminal is accordingly 
adjusted. To implement this timing method, the ephemeris data should be available at a 
reference terminal and every traffic terminal. Since the accuracy of the generated TDMA 
timing depends on the accuracy of the ephemeris data, the ephemeris data needs to be 
periodically updated to maintain a sufficient timing accuracy. Furthermore, when a terminal 

20 loses a reference burst, a separate transmission link should be established to update the 

ephemeris. Associated data processing and fransmission for this updating increase the system 
complexity and operational cost. 

In a TDMA satellite communications system, the timing of every traffic terminal is 
typically adjusted with reference to the timing of a reference terminal, as shovra by way of 

25 example in FIGURE 4. During normal operation, the receive timing of each traffic terminal 
is derived from reception of fraffic bursts transmitted from the reference terminal. The 
transmit synchronization is typically achieved by a feedback control process which includes 
measuring timing offset of a traffic burst at the reference terminal, sending the timing 
correction infoimation to the traffic terminal, and adjusting the burst transmit timing at the 
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traffic terminal. A reference burst is used for sending the correction information. If, for any 
reason, the traffic terminal loses a reference burst, neither the receive timing nor the 
correction information for adjusting the transmit timing is available at the traffic terminal. 
Thus, there is a need for a method that can be used to support accurate continued system 
5 operation by allowing continued generation of the receive and transmit TDMA timings at a 
traffic terminal, or the receive TDMA timing at a reference terminal, or when the receive 
reference burst is lost. 

SUMMARY OF THE INVENTION 

In view of the deficiencies of conventional systems, it is an object of the present 
1 0 invention to maintain accurate generation of transmit and receive timings at a traffic terminal 
^ regardless of whether a receive reference signal is lost. It is an additional object of the present 

invention to maintain accurate timing at a traffic terminal without a need for an external 
inputting of satellite ephemeris data. 

The present method and apparatus are appUcable to conventional user traffic terminals 
15 that normally receive their timing synchronization by accessing a reference terminal to 
correct their individual flj^heel timings, and to a flywheel timing generation within the 
reference terminal. A flywheel timing for satellite communications synchronizes TDMA 
communications. The flywheel momentum maintains a given accuracy for a period of time 
with pulses based on the overall rate of reception. The present invention allows the continued 
20 generation of the receive and transmit TDMA timings at a traffic terminal, or the receive 
TDMA timing at a reference terminal, or when the receive reference burst is lost. After an 
initialization period of normal operation, the present method does not require any external 
information such as satellite ephemeris, and can maintain sufficient timing accuracy for 
; several hours of flywheel operation. Delay values associated with generating the timings are 

25 predicted based on prior drift history, which is generated based on measurements at the 

terminal diuing normal system operation. The delay values can be calculated according to a 
number of different formulas, and various circuit implementations are envisaged. For 
example, in an INTELSAT system such as that disclosed in IEEE Journal on Selected Areas 
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in Communication . Vol. SAC-1, No. 1, pp. 165-173, an accuracy of at least +/- 14 symbols is 
achieved when measuring over a two-hour flywheel operation. 

The present flywheel timing generation method uses the prior history of TDMA 
timing drift measured at the terminal during normal operation. 
5 The reference terminal derives the transmit timing from its local timing source. 

However, the reference terminal derives the receive timing from reception of reference 
bursts. The present method can be used to generate the receive flywheel timing when the 
reference terminal loses the reference burst. 

The satellite range change due to the orbit inclination and eccentricity is periodic with 

10 one sidereal day. If an inclination and eccentricity do not change, a flywheel delay value can 
be obtained by a normal delay value measured one sidereal day prior to the flywheel 
operation. However, gravitational pull of the sun and the moon at times can change the orbit 
inclination by 0.005% per day. Unless there is compensation, the range change due to the 
inclination change may cause unacceptable timing offset in TDMA system operation. For 

15 example, the maximum range difference reaches approximately 420 meters, or 170 symbol 
periods in the INTELSAT TDMA system, between transmit and receive earth stations, 
assuming that the earth stations are located at 45*^ latitude and 45** longimde from the satellite 
nadir direction. Therefore, to achieve sufficient accuracy, the present method accounts for 
periodic satellite range change due to the orbit inclination and eccentricity as well as the daily 

20 change of the orbit inclination. 

The present method includes measuring a satellite drift in the north/south direction at 
an earth station, generating a history of the measured drift over a period of time, and 
generating flywheel timing values based on delay values predicted according to the measured 
satellite drift history over a predetermined time. The predetermined period of time can be one 

25 sidereal day. The history of the measured drift can be maintained at the reference terminal as 
well as at the traffic terminals, so that flywheel timing values can be generated either 
independently or by sharing of information. The flywheel timing values, for a TDMA 
sateUite communications system having a traffic terminal and a reference terminal, are 
generated by calculating a number of symbols with respect to reference pulse timing of the 
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traffic terminal. During normal operation, the reference terminal transmits correction 
information to the traffic terminal. The correction information can include a reference burst 
used by the traffic terminal to derive a receive timing. Synchronization for transmitting can 
be achieved by measuring timing offset of a traffic burst at the reference terminal, and 
5 adjusting a burst transmit timing at the traffic terminal. 

An apparatus according to the invention can include a computer system having a 
processor and a memory, the memory including software instructions adapted to enable the 
computer system to perform the steps of: generating a reference pulse stream with a period of 
one control firame; measuring and recording a plurality of receive time delay and transmit 

10 time delay values for a satellite communication signal over a predetermined period of time; 
designating, for every control frame interval, start of receive fi-ame delay and start of transmit 
firame delay values based on the control firame period and based on the recorded time delay 
values, referenced to a designated time, from a designated value for receive frame delay; 
generating a flywheel receive start timing by counting a calculated number of symbols fi-om a 

1 5 corresponding designated reference pulse; and, from a designated value for start of transmit 
fi-ame delay, generating a flywheel transmit start timing by counting a calculated number of 
symbols fi-om a corresponding designated reference pulse. 

The computer system can include any means for generating flywheel timing values 
considering either the daily inclination change, or a daily delay value change, due to a 

20 satellite drift in the north / south direction. The daily delay value is computed as a fimction of 
a maximum time difference due to the satellite drift in one sidereal day. 

A circuit for flywheel operation in a satellite communications system, according to 
the present invention, includes a first counter that measures a receive delay time during a 
normal operation, the first counter operative to receive a predicted receive delay value and 

25 generate a flywheel receive control timing during flywheeling operation, a second counter 
that measures a transmit delay time during a normal operation, the second counter operative 
to receive a predicted transmit delay value and generate a flywheel transmit contiol timing 
during flywheeling operation, a first latch operative to record the measured receive delay 
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time, and a second latch operative to record the measured transmit delay time. The circuit can 
have a symbol clock operative to generate a reference pulse stream. 

BRIEF DESCRIPTION OF THE DRAWING 

5 In the accompanying drawing: 

FIGURE 1 illustrates the TDMA timing during system operation according to an 
embodiment of the present invention; 

FIGURE 2 illustrates TDMA control frame designation for flywheel operation 
including designation of Start of Receive Control Frames (SORCF) upon loss of receive 
10 timing, according to an embodiment of the invention; 

FIGURE 3A illustrates a circuit configuration during normal operation; 
FIGURE 3B illustrates a flywheel circuit implementation during flywheeling 
operation; 

FIGURE 4 illustrates a general architecture for a satelhte communication system that 
15 includes a reference terminal sending periodic timing correction information to traffic 
terminals; and 

FIGURE 5A illustrates an initialization used in an embodiment of a method 
according to the present invention; 

FIGURE 5B illustrates a normal operation phase used in an embodiment of a method 
20 according to the present invention 

FIGURE 6 illustrates an operation upon loss of receive timing, used in an 
embodiment of a method according to the present invention; and 

FIGURE 7 illustrates a terminal apparatus according to an embodiment of the present 
invention. 

25 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
FIGURE 1 illustrates the TDMA timing during system operation. The receive timing 
seen at a traffic terminal drifts over time as the range between the satellite and traffic terminal 
changes due to satelUte motion. In this case, a "control fi-ame" is defined as an integer 
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multiple of a TDMA frame. Timing is measured and adjusted on a control frame basis. 
During normal operation, receive timing is generated by receiving the reference burst, and 
transmit timing is generated using a delay value which is provided by the reference terminal. 
For example, the transmit timing at the i*^ control frame is generated using Dp[i], a delay 
5 value provided by the reference terminal for the i* control frame. If the traffic terminal loses a 
reference burst, the present method generates i*^ confrol frame timing using an SORCF (Start 
of Receive Control Frame) delay value DrP] and an SOTCF (Start of Transmit Control 
Frame) delay value DjCi]. The Reference Pulse Timing is represented by RPT. The delay 
values are generated using various formulas. The delay value DR[i] represents the number of 

10 symbols with respect to the reference pulse tuning of a fraffic tenninal at the i* control frame. 
DR[i] and Dj{i] change due to the range change and due to the drift of the traffic tenninal 
clock with respect to the reference terminal clock. 

Using the symbol clock, a reference counter at the traffic terminal generates a 
reference pulse stream with a period of one control frame. The approximate location of the 

15 reference pulse is estabhshed relative to the SORCF to facilitate the fl5wheel timing 

generation. An exemplary embodiment of a flywheel timing generation method is shown in 
FIGURES 5 A, 5B-, and 6. At network initialization, the reference pulse is appropriately 
placed such that Dr^] and DT[i] values do not exceed the number of symbols in one confrol 
frame during the system operation. Step 101. During normal operation, DR[i] and DT[i] are 

20 measured and recorded, Step 201. A time stamp is put on each DrP] and DT[i] value. Step 
202. A file of DR[i] values is created for one sidereal day (86164.091 seconds) plus 4NtTc 
seconds, Step 203, where Tc denotes a confrol frame period in seconds and denotes an 
integer parameter that is optimized for a required flywheel duration. The integer can be 
experimentally derived and changed either by an operator or automatically according to a 

25 history file, or according to a predetermined condition. Typical values for 4NxTc are around 
7200 seconds and do not exceed 21600 seconds. A file of DrCi] values is created for one 
sidereal day, Step 203. In the event of loss of receive timing, SORCF delay values are 
designated for fljnvheel operation. FIGURE 1 and FIGURE 2 illustrate an example of this 
designation of SORCF delay values. A measured SORCF delay value Dr[0], obtained just 
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prior to the timing loss, is used as a reference, Step 301. From this position Dr[0], previous 
measured values and subsequent flywheel values are designated as DrH, (i = -4N,, . . . -2, - 
1, 0, 1,2,.. .) in every control frame interval, Step 302. A measured SORCF value, obtained 
one sidereal day prior to the mstant of the Dr[0] measurement, is designated as Dri[0], Step 
5 303. From this position at Dri [0], previous and subsequent measured values are designated as 
DRi[i], (i = -4Nx, . . . -2, -1, 0, 1, 2, . . .) in every control frame interval, Step 306. 
A measured SORCF delay value, obtcdned just prior to the timing loss, is designated as 
Dt[0], Step 304. From the position at Dt[0], subsequent flywheel values are designated as 
I^T[i], (i = -4Nx, ... -2, -1, 0, 1, 2, . . .) in every control frame interval. Step 306. A measured 

10 value, obtained one sidereal day prior to the instant of the Dt[0] measurement, is designated 
as Dti[0], Step 305. From this position at Dti[0], previous and subsequent measured values 
are designated as DTi[i], (i = -4Nt, ... -2, -1, 0, 1,2,.. .) iii every control frame interval, Step 
306. The flywheel values used to generate receive timing are calculated using a flywheel 
timing generating equation derived on the basis of a model that approximates the daily delay 

15 value change due to the sateUite drift and that assumes a satellite range change being 

sinusoidal with one sidereal period. Several non-limiting variations for a flywheel timing 
generating equation are detailed below. It is understood that one skilled in the art can use 
approximations other than those illusfrated for deriving alternative computations for a daily 
delay value change. In a first example, using the recorded measured values DRi[i] and DR[i], 

20 the flywheel values DR[i] (i = 1, 2, . . .) are calculated, Step 307, by 

DR[i] = DRi[i] + Dr[-2NJ - Dri[-2NJ 
+ INT[(Tsd/ 471 N.To) { Dr[-N.] - Dri[-N.] - Dr[3NJ + Dri[3NJ} sin{27r (Tc/Tsd)(i + 2N,)}] 
+ INT[(Tsd/ An N,T,)' { Dr[0] - Dri[0] - Dr[-4N,] + Dri[-4N,] - 2Dr[-3N,] + 2Dri[-3N,]} 
25 .[1 - cos{2:t (Tc/Tsd)(i + 2N,)}]]. 

where Tsd, Tc, and INT["] respectively denote one sidereal day duration, confrol frame period, and 
the nearest integer. For example, rNT[1.6]=2; INT[1.3]=1; and rNT[1.5]=2. Using the recorded 
measured values DTi[i] and DT[i], the flywheel values DT[i] (i = 1, 2, . . .) can be calculated. Step 
30 308, for example, by 
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DtP] = Dnti] + Dt[-2N,] - Dti[-2NJ 
+ INT[(Tsd/ 471 NtTc) { Dt[-NJ - Dt][-N,] - Dt[3Nx] + Dti[3NJ} sin{27i (Tc/Tsd)(i + 2N,)}] 
+ INT[(Tsd / 471 N.Tcf { Dt[0] - Dti[0] - Dt[-4NJ + Dti[-4Nx] - 2Dt[-3NJ + 2Dti[-3N,]} 

.[l-cos{27i(T/r,d)(i + 2N0}]]. 

5 A derivation of these exemplary equations is shown below. From the obtained values 

for DrP] and DT[i], the traffic terminal generates receive and transmit timings. Start of i* 
SORCF-F (flywheel SORCF) is generated by counting DR[i] symbols from the i* reference 
pulse, Step 309. Start of i*^ SOTCF-F (flywheel SOTCF) is generated by counting DT[i] 
symbols from the i'*' reference pulse, Step 310, 

1 0 It will understood by one of ordinary skill in the art that the exemplary method 

described herein can be implemented in part or as a whole, and that the identification of 
individual steps is non-limiting. A use of the invention in any appropriate programming 
language is envisaged, and the method may be employed by any suitable hardware, including 
a stand-alone terminal or processor, multiple processor configuration, or in a network. 

15 A terminal utilizing the present method can be a simple telephone terminal or a 

complex multichannel system including local area networks (LAN) or wide area networks 
(WAN) within its control. 

FIGURE 7 illustrates a terminal 1 used to implement a flywheel timing generation 
method. The terminal 1 can include a computer system having a processor 2 and a memory 3, 

20 the memory 3 including software instructions adapted to enable the computer system to 

perform the steps of: generating a reference pulse stream with a period of one control fi-ame; 
measuring and recording a plurality of receive time delay and transmit time delay values for a 
satellite communication signal over a predetermined period of time; designating, for every 
control frame interval, start of receive frame delay and start of transmit frame delay values 

25 based on the confrol frame period and based on the recorded time delay values, referenced to 
a designated time, from a designated value for receive frame delay; generating a flywheel 
receive start timing by counting a calculated number of symbols from a corresponding 
designated reference pulse; and, from a designated value for start of transmit frame delay, 
generating a flywheel transmit start timing by counting a calculated mmiber of symbols from 

30 a corresponding designated reference pulse. The time delay measurement circuit 4 can 
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comprise counters and latches, as shown in FIGURES 3A, 3B. The pulse generator 5 can 
supply a system clock, a reference pulse stream, flywheel receive start timing pulses, 
flywheel transmit start timing pulses, synchronization words, as well as timing control 
information. 

5 The computer system will tj^ically have a 32-bit, or larger, microprocessor interface, 

with a data interface able to accept corrrniunications that include serial, parallel, and 
synchronous or having a timecode interface. The computer system also has a data output 
adaptable to generating j&ame data structures that include flywheel frames. The computer 
may be implemented using very large scale integrated circuit (VLSI) components of various 

10 protocols, and can be arranged as subsystems that are specific to a certain data format or 
optimized for various error detection, coding schemes, or transfer rates. 
Flywheel Circuit Implementation 

An exemplary embodiment of a circuit implementation for flywheel operation is 
illustrated in FIGURES 3A and 3B. The circuit operation is described as follows. 

15 Using a symbol clock of a traffic termmal, a reference pulse stream with a period of 

one control frame is generated. During normal operation, the circuit is configured as shown in 
FIGURE 3 A. Dr and Dt values are measured by the respective UP counters 11, 12. Then, 
using the respective latches 21, 22, Dr and Dt values are recorded. During flywheel 
operation, the circuit is configured as shown in FIGURE 3B. The predicted Dr and Dt 

20 values are sent to the respective DOWN coimters 31, 32. Then, from these DOWN counters, 
the SORCF-F and SOTCF-F flywheel timings are generated. The symbol Nc represents the 
number of symbols in one control frame. 
Flywheel Timing Generation Example 

To better illusfrate the present method, an example is now provided for generating 

25 flywheel timing in an application to the INTELSAT TDMA system. For two hours of 
flywheel operation, a simulation indicates that an optimum value of N, is equal to 1800. 
Substituting N^ = 1800 into the above equation for determining the flywheel values DR[i], and 
rearranging terms provides an SORCF fl)wheel generating equation: 
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DrW = DriH + Ai 4- INT[A2 sin{7.467 x 10"^ (i + 3600)} 
+ A3(l - cos {7.467 X 10"^ (i + 3600)})], 

where 

Ai=Dr[0]-Dri[0], 
5 A2 = 3.72(Dr[-1800] - Dri[-1800] - Dr[-5400] + Dri[-5400]), and 

A3 = 13.8385(Dr[0] - Dri[0] + Dr[-7200] - Dr,[-7200] - 2Dr[-3600] + 2Dri[-3600]); 

likewise, an SOTCF fl3nvheel value generating equation becomes: 
DtH = DtiW + Bi + INT[B2 sin {7.467 x 10"^ (i + 3600)} 
+ 83(1 - cos {7.467 X 10"* (i + 3600)})], 

10 where 

Bi=Dt[0]-Dti[0], 

B2 = 3.72(Dt[-1800] - Dti[-1800] - Dt[-5400] + Dti[-5400]), and 

B3 = 13.8385(Dt[0] - Dti[0] + Dt[-7200] - Dti[-7200] - 2Dt[-3600] + 2Dti[-3600]). 

15 In such an example, the inventors' experimental computer simulation results, which 

include quantization and approximation errors, indicate that an acciu-acy of +/- 14 symbols is 
achieved during a two-hour flywheel period for both SORCF and SOTCF. 
Derivation of a Flywheel Value Generating Equation 

Referring to FIGURE 1, a SORCF delay value changes in accordance with the range 

20 of the reference terminal to satellite plus the satellite to traffic terminal. A satellite range 
change due to orbit inclination and eccentricity is periodic over time with one sidereal day 
period. Therefore, if the orbit inclination and eccentricity do not change in a day, a delay 
value at a certain time can be given by that of one sidereal day prior to that time. However, 
the orbit inclination changes during a day. Therefore, the flywheel value generating equation 

25 is derived considering the daily inclination change due to the satellite drift in the north/south 

direction. The east/west satellite drift is ignored because the impact on the range change is 

much smaller than that of the north/south drift. 

Since a satellite range change is sinusoidal with one sidereal period, the daily delay 

value change due to the satellite drift, Z(t), can be approximately modeled by: 

30 Z(t) s DrO) - DRi(t) 

= Ki sin(27tt / Tsd + 6) + K2, 

where: 

DR(t) denotes the time difference between the reference pulse and SORCF, 
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DRi(t) denotes the time difference between the reference pulse and SORCF at one sidereal 
day prior to a time t, 

Ki denotes the maximum time difference due to the satellite drift in one sidereal day, 
K2 denotes the timing difference due to the traffic terminal clock drift with respect to the 
5 reference terminal clock, 

Tsd denotes one sidereal day, i.e., 86164.091 seconds, and 

9 denotes a random phase associated with the daily delay change. 

Using this definition of the daily delay value change due to the satellite drift, Z(t), a 
predicted value, DR(to + At), can be written as: 
10 DR(to + At) = DRi(to + At) + DR(to) - DR,(to) + Z(to + At) - Z(to). 

The tenn Z(t + At) - Z(t) can be written as: 

Z(t + At) - Z(t) = (Tsd / 2n) sin(27iAt / T^d) d/dt Z(t) 
- (Tsd / 27t)^{1 - cos(27iAt / Tsd)} dW Z(t). 
15 For At « Tsd, 

d/dt Z(to) « {Z(to + At/2) - Z(to - At/2)} / At, 

and 

dW Z(to) {d/dt Z(to + At/2) - d/dt Z(to - At/2)} / At 

« {Z(to + At) + Z(to - At) - 2Z(to - At/2)} / AtI 

20 

In addition, Z(to + At) - Z(to - At) can be represented as follows: 

Z(to + At) - Z(to - At) « (Tsd / 27cAt) {Z(to + Ax/2) - Z(to - At/2)} sin(27tAt / T^d) 
+ (Tsd / 27cAT)^{Z(to + At) + Z(to - At) - 2Z(to - At/2)} (1 - cos(27tAt / Tsd)} • 

25 Putting the representation for Z(to + At) - Z(to - At) into the predicted value DR(to + At) and 
using the definition of Z(t), the flywheel timing generating equation is given by: 



DR(to + At) = DRi(to + At) + DR(to) - DR,(to) + (T^d / 27rAT) {Z(to + At/2) 

- Z(to - AT/2)}sin(27rAt / Tsd) + (Tsd / 27cAt)^ {Z(to + At) + Z(to - At) 
30 - 2Z(to - At/2)} { 1 - cos(27iAt / Tsd)} 

= DRi(to + At) + DR(to) - DRi(to) + (Tsd / 27cAt) {DR(to + At/2) 

- DRi(to + At/2) - DR(to - At/2) + DRi(to - AT/2)}sin(27cAt / T^d) 

+ (Tsd / 27tAT)^{DR(to + At) - DRi(to + At) - DR(to - At) + DR,(to - At) 

- 2DR(to - At/2) + 2DRi(to - At/2)} {1 - cos(27rAt / Tsd)}. 

35 
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By using this equation, DR(to + At) can tiius be calculated from DR(t), t < to + At, and DRi(t), t 
< to + Ax. In other words, DrO), t < to + At, can be predicted using previous values of DR(t) 
and DRi(t), assuming that At is not greater than Tsd + Ax. At a traffic terminal, the time 
difference between the reference pulse and SORCF is measured once every control frame. 
5 The time difference can be expressed, for example, as a nimiber of symbol periods. In 

association with the flywheel timing generating equation, the measured values are defined by: 

DR[i] s INT[DR(to + At + i-Tc) / T, ], 
DRi[i] = INT[DRi(to + At + i.Tc) / Tj ], and 
10 At s 2.Nt«Tc, 

where Ts is the symbol period. Then, the SORCF flywheel timing generating equation can be 

written as: 

DR[i] = DriH + Dr[-2NJ - Dri[-2N,] + INT[(Tsd / 27iN.Tc){Dr[-N,] - Dri[-NJ 

- Dr[3Nx] + Dr,[3NJ} sin{27i (J/Tsd)(i + N,)}] + INT[(Tsd / 2% N.T^f {Dr[0] - Dri[0] 
15 - Dr[-4Nx] + Dri[-4NJ - 2Dr[-3N,] + 2Dri[-3N,]} [1- cos {2n (Tc/Tsd)(i + N,)}]]. 

This algorithm compensates flywheel timing offset due to the traffic terminal clock drift with 
respect to the reference terminal clock since the term DRi(to + At) - DRi(to) in the flywheel timing 
generating equation compensates the time offset generated for a duration of At. 

20 The SOTCF delay value can also be accurately deduced. For TDMA synchronization, the 

timing seen at the satellite should be synchronized. The SOTCF delay value changes in accordance 
with the range of the reference terminal to satellite minus the range from the satellite to traffic 
terminal. Since the range change is sinusoidal over time with one sidereal day period, the daily 
SOTCF delay value due to the satellite drift in the north / south direction can be determined by using 

25 the same daily delay value change due to the satellite drift, Z(t), modeled as noted above. Therefore, 
by using the SORCF flywheel timing generating equation just described for the SORCF delay value, 
the SOTCF flywheel value generating equation can be written as: 

DxEi] = DTi[i] + Dt[-2NJ - Dti[-2NJ + INT[(Tsd / IttN.Tc) {Dt[-N,] - Dti[-N,] 
30 -Dt[3NJ + Dti[3N,]} sin{27:(Tc/Tsd)(i+N.)}] + INT[(Tsd/27rN,Tc)' {Dt[0] - Dti[0] 

- Dt[-4N,] + Dti[-4N,] - 2Dt[-3N,] + 2Dti[-3N,]} [1- cos {In (To/Tsd)(i + N,)}]]. 
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Thus, a flywheel timing generation method, circuit, and computer system are provided 
for a satellite communications system. One skilled in the art will appreciate that the present 
invention can be practiced by other than the described embodiments, which are presented for 
purposes of illustration and not limitation, and the present invention is limited only by the 
claims that follow. 
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